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A dynamical model for density fluctuations in a one-component fluid near the critical point is pro- 
posed and compared with existing measurements in carbon dioxide. The model is described by a set 
of linearized hydrodynamic equations modified to include a nonlocal pressure density relationship and 
to include relaxation in the volume viscosity. Parameters for the model are found which are consistent 
with bulk measurements. With these parameters the model reproduces, within experimental un- 
certainty, the observed Brillouin spectrum of critical opalescence in C0 2 . The low frequency volume 
viscosity is found to diverge as the— 1/3 power of T— T r . An additional modification of the hydrodynamic 
equations, a frequency dependent thermal conductivity, is considered, but no definite conclusions can 
be reached as the Brillouin line-width data lack sufficient precision. 
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1. Introduction 

The dynamical structure of density fluctuations in a 
fluid is directly related to the spectral structure of 
light scattered by the density fluctuations. This struc- 
ture has been the subject of a large number of theo- 
retical [l] 1 and experimental [2~4| papers. It has been 
established that the linearized hydrodynamic equations, 
suitably modified to include internal degrees of 
freedom, represent a satisfactory model for the 
dynamics of the long wave length density fluctuations 
in many simple fluids. 

In this paper we propose to compare a model for the 
dynamics of density fluctuations in a one-component 
fluid near the critical point with recent experimental 
observations of the spectral structure of critical 
opalescence. The intent is to determine the main 
features of a model for critical fluctuations. We are 
not so much concerned with the precise values of 
parameters occuring in the model as we are that these 
parameters be in agreement with the known properties 
of the fluid system. The measurements are for carbon 
dioxide at the critical density for temperatures greater 
than the critical temperature. Three modifications of 
the linearized hydrodynamic equations will be ex- 
amined. The first modification is the introduction of 
a nonlocal relationship between density and pressure 
fluctuations [5]. The second is the introduction of a 
frequency dependent volume viscosity and the third 



modification is the introduction of a frequency depend- 
ent thermal conductivity [6, 7]. For convenience we 
assume the single relaxation time form for the fre- 
quency dependence of the volume viscosity and the 
thermal conductivity. Given the model and the experi- 
mental results currently available we conclude that 
the first and second modifications in the hydrodynamic 
equations are required. It is not possible to decide 
whether or not the thermal conductivity is frequency 
dependent. 

2. Dynamical Model 

Before we introduce the model for the density 
fluctuations in terms of a set of linearized hydro- 
dynamic equations let us briefly state what is found 
experimentally. The spectrum consists of an intense 
central component and two relatively weak Brillouin 
components. The width of the central component may 
be expressed (in units of rad/sec) as 



r c = x k-[i+ek 2 ] 



(1) 
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where x = VpoCp is the thermal diffusivity at constant 
pressure. Here X is the thermal conductivity, C p is the 
specific heat at constant pressure and po is the number 
density. The wave vector of the fluctuation is k. For the 
experiments to be considered /r=2.18Xl0 7 m 1 . 
The temperature dependence of \ is found to be, for 
CO2 above the critical point, proportional to 
I T—T c I ~° ■ 73±0 - 02 where T c is the critical temperature 
(T c — 303.2 K for C0 2 ). Swinney and Cummins found 
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£=(2.4±0.51)Xl0- 9 /|r-r c | 2 / 3 m for C0 2 at the 
critical density [4]. 

The Brillouin components exhibit a large dispersion 
in the frequency shift o)B = VBk; that is, vb>Co the 
adiabatic sound speed. At the same time the width 
of the Brillouin components is appreciably greater 
than the width predicted by the unmodified hydro- 
dynamic equations [2, 3]. 

Finally, the ratio of the intensity of the central com- 
ponent to the intensity of the Brillouin components is 
found to be proportional to \T-T C \ -^±0.03 [2 , 3]. 

The dynamical model for density fluctuations of 
wave vector k is 



dp k (t) 
dt 



-po<MO=0 



(2) 



W>M = CK (1 + w ) {t) +CMr k 2 Tk{t) 

dt poy y 

-b*k*Mt)-k 2 f'bWfait-t'W, (3) 

(5t dt dt 



k> {' A(t' 
Jo 



)T k {t-t')dt'. (4) 



Here pk(t) is the A:th spatial Fourier component of 
the density given that at t = that component had the 
value ptc.Tk(t) and ift k (t) = ik • v k (t) are the cor- 
responding Fourier components of the temperature 
and velocity field. In eq (4) the entropy has been re- 
placed by the temperature and the density using local 
thermodynamic equilibrium arguments. The pressure 
term in eq (3) has been also replaced by the tempera- 
ture and density using a nonlocal relationship between 
the pressure and the density. The Fixman modifica- 
tion [5] was used so that 

p»(0-= 3 (i+f/K-)P»(0+=^ e r»w. 



Here k~ 1 is the two particle correlation length, po is 
the equilibrium number density and m is the mass of 
a molecule. Other quantities are Co, the low frequency, 
adiabatic sound speed, y = C P IC v is the ratio of the 
specific heat at constant pressure C p to the specific 
heat at constant volume C v and ft T is the isobaric 
thermal expansion coefficient. The nonrelaxing longi- 
tudinal kinematic viscosity is bo= (4>/3r) s -\-r) v )lmpo 
where r) s is the shear viscosity and r) v is the frequency 
independent volume viscosity. The relaxing part of 
the volume viscosity is contained in the convolution. 
We assume 

b(t) = (Cl-Q)exp(-tlr). (5) 

The time dependence of b is arbitrarily chosen to be 



an exponential decay. While there is no reason to 
expect this simple form, we shall see that the data 
are, at present, adequately represented using this 
form. The "high frequency" speed of sound is C x . 
Equation (5) is a frequently used form for relaxing 
viscosities. 

In eq (4) \ is the static thermal conductivity and 
the relaxing part is A(t). Again we choose an expo- 
nential decay: 



rA (t) = (Xoo- k ) exp (- tfr) . 



(6) 



The high frequency thermal conductivity X^ is prob- 
ably equal to the thermal conductivity when T—T c 
> 20 K. Except near the critical point the thermal 
conductivity is insensitive to small changes in the 
temperature [8] . The form for eq (6) was chosen on the 
supposition that only the critical anomaly in the ther- 
mal conductivity can exhibit relaxation. Again the 
factor exp [— tfr] was chosen for simplicity. 

The spectrum of light scattered by density fluctua- 
tions is specified by the correlation function 
( p k (ct))p- k ) where p k (a)) is the Fourier transform 
of p k (t). We shall determine (p k (co)p- k ) for the 
model by solving eqs (2)-(6) as an initial value problem 
for p k (t) given that p k (0) = p k . Put another way, we 
use the equations of motion (2)-(6) in lieu of the 
conditional probability which relates p k (t) to p k . If 
desired, the ensemble average indicated by the angular 
brackets (...) may be computed using an equi- 
librium ensemble. The initial value (p k p- k ) is not of 
special interest here as the Fourier components are 
independent in this model. [If there were coupling 
between p k and p k+q the initial values would be im- 
portant.] The reader is referred to reference [1] for 
more detailed discussion of these calculations. 

Solution of the model is straightforward when the 
Laplace transform is employed. The transform of 
pk(t) is 



Pk(z) 



-I' 



e-"p k (t)dt 



(7) 



with similar expressions for T k (z) and ijj k (z). With 
the assumption that ( p k T- k ) = ( p k ^- k )=0 it 
follows that 



(pk(z)p- k ) = F(z) 
<|p*l 2 > G(z) 



(8) 



with 



F(z)=z* 



+z 



1 + 



(a x —a )k-T 
1 + zf 

(Q-CfikH 

l+ZT 



+ 



b k 2 + a k 2 + 
| (a,-ao)Af / | (Cl-Cj)kH \ 

1+ZT \ 1+ZT / 

(l--)c 2 k 2 + a k 2 (b k 2 + 



(9) 



(Cj-Q)k 2 r 

l + ZT 
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and 



G(z)=z>h+ {a ~r°°}^] 

L 1+ZT J 



+ Aa»k 2 +b () k 2 + 
(a x — a )k 2 r 



+ - 



1+zr 



(C 2 X -C§)A* T 

1 + ZT 

(C£-CJ)**r 



W 2 +" 



1-hZT 



L 1+ZT 

+ C 2 A: 2 (l + A: 2 /r>) 

+ ^(1 + P//C 2 ) (^- Q o)^ 2 ^ 



1+ZT 






(10) 



We have introduced the quantities a = A /poC r and 

An exact expression for (p k (co)p- k ) may be ob- 
tained from eq (8) by replacing z by io) and taking 
1/7T times the real part of the resulting expression. 



_ <p k ((x))p- k > F(ico) 

77 — ; — — ne -ftt. — v-. 

When this is done we find 

< pkMp-k > F X G X +F>G 2 



<Ipa-| 2 > (C,) 2 + (C 2 ) 2 

whereF(i'w) =F\ + iF 2 andG(io)) — G\ +1G2. 
The explicit forms for F\ etc. are 



(11) 



(12) 



Fi(oo) = — (o 2 
(ayo — ao)k 2 f I 

1 + (x) 2 f 2 \ 



1 + 



(a x -ao)k 2 f {Cl-Cl)k 2 r 2 



l+o> 2 f 2 1 + co 2 t 2 

(a^ — ap)k 2 f ( {Cj—Co)k 2 T(r-\-f) 



1+ W 2 T 2 



-+b k 2 f 



1 -±) CgJfc 2 + a Jfc 2 (m 2 + (C \ + ^ 2T > 



r( x ^ (o»- ao)k 2 r 2 [ 1 
F 2 (co)=g> 3 — „_„ — 1 

1 +(ti*T- 



(Cj-C' 2 )k 2 r 2 

1 + W 2 T 2 
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(a x — a { ))k 2 f 
l + co 2 f 2 

— a k 



L J 2 1 7 2 _!_ (^x ^o)k"T 
bok z + aok~ + - — 



1 + 0) 2 T 2 



\ 1 + arr- 



„ (C*-Cg)AV 



1 + wV 2 



(14) 



Gi(w) = — a> 4 



l+w 2 f 2 V l + to 2 T 2 



(aoo — a )A: 2 f ^ l9 , (C 2 -Cg)A: 2 T 
l + w 2 f 2 



M 2 +- 



: 2 -Cq 2 )A; 2 t \ 
1 + a> 2 r 2 / 



,, (C 2 -C 2 )Fr 2 C 2 k 2 L.k 2 
■ ao/c z — t - : — r^ — 



l + wV 2 



y 



(aoo — a )k 2 T 2 



l + w 2 f 2 



+— cik 2 {\ + k 2 iK 2 ) 
y 



(15) 



and 

Go(o)) = — CO 3 



1 + 



(a x -a )tff) -(C 2 -C 2 )/c 2 t 2 



l + co 2 r 2 



1 + W 2 T 2 



(a x -aoU 2 r A rr| (C'-CjQ/tMt+t ) 



C 2 )Fr(T+f )\ 
l+co 2 r 2 VT U " ' ' 1 + 0) 2 T 2 /. 

y V /< 2 / l + co 2 r 2 J* 



(16) 



3. Comparison With Experiment 



The frequency spectrum of density fluctuations 
for the model may be computed by evaluating eq (12) 
for various values of a>. To do this requires specified 
values for the parameters ao, a x , bo, Co, C x , y, k, r, 
r, and A. The wave vector k is determined by the 
experimental configuration. In what follows we shall 
use A = 2.18 X 10 7 m~\ the value appropriate to the 
measurements of Ford, Langley, and Puglielli [3]. 
It is possible to make reasonable estimates for a x , 
60, Co, and C x using various experimental results. 
We can also specify a /y reliably. The remaining 
parameters k, t, f and y are to be determined by 
requiring (p k (oj)p-k) to correspond to the spectrum 
observed for C0 2 . 

First let us specify a x and 60. The critical anamoly 
in the thermal conductivity may be detected 40 K 
above the critical point in C0 2 . However, most of the 
increase occurs for \T— T c \ < 20 K so a x was chosen 
to be 1 X 10~ 7 m 2 /s which is comparable to a at T—T ( 
= 20 K [8]. The shear viscosity and the nonrelaxing 
part of the volume viscosity are contained in 60. 
A value 1.5 X 10 _7 ra 2 /s was chosen by assuming 
r] r ~ rj s and using measured values for tj s [9]. 

The low frequency sound speed values Co were taken 
from figure 2 of reference [2]. These values correspond 
to the thermodynamic sound speed corrected for the 
vibrational degrees of freedom which are "frozen out" 
at the Brillouin frequency. This is the only place where 
the vibrational degrees of freedom are expected to be 
significant. Near the critical point the configurational 
contribution to the specific heat is large compared to 
the vibrational contribution so y, the ratio of the 
specific heats, is not expected to depend significantly 
on the vibrational specific heat. The high frequency 



362-197 0-69— 3 



595 



sound speed values C x were assumed to be slightly 
larger than the speeds obtained from figure 3 of refer- 
ence [3]. For T—T c < 1 K, C^ = 194 m/s was found to 
yield good agreement with experiment. Values of C x 
for T— T c > 1 K are shown in table 1. 



Table 1 



T-T c , K 


Co, m/s 


C x , m/s 


0.05 


138 


194 


.1 


140 


194 


1.0 


160 


194 


2.0 


177 


196 


5.0 


198 


208 


10.0 


225 


225 



The low frequency sound speeds, Co, were taken from figure 2 
of reference [2]. The high frequency sound speeds, C x , were taken 
from figure 3 of reference [3J. 



The other quantity which may be specified is 
ao/y = \{)lp()C p the thermal diffusivity at constant 
pressure. The width of the central component is, to 
lowest order in A; 2 , given by kok 2 lpoC p (see eq (1)). The 
k 4 terms of this model may be found by seeking small 
solutions to the dispersion equation 



G(z)=0; 



(17) 



G(z) is given in eq (10). The small z solution to eq(17)is 
aok 2 



y L 



i + 



F/P-ao^r(-p-l)-^^]. (18) 



Terms which are of order 1/y compared to the retained 
terms have been dropped. Comparison with eq (1) 
indicates that 



^- aoT (§ _1 



y 



(19) 



Swinney and Cummins have determined a /y by 
measuring the width of the central component in 
C0 2 and extracting the k 2 coefficient from those 
measurements [4]. 

The remaining parameters were evaluated as follows. 
An initial value was assumed for y. [Once y is speci- 
fied, a () is obtained from the data of ref. 4]. The vis- 
cosity relaxation time r was adjusted so that the 
Brillouin line width T H , as determined from the spec- 
trum obtained by evaluating eq (12), was in agreement 
with experiment. The sound speed v{k) is not very 
sensitive to the value of r as (jl)bt = v (k) kr > 3. Next 
the values of k and f were adjusted so that the central 
component line width is reasonable. In practice 
"reasonable" means £ 2 nonnegative. The intensity 
ratio was then checked for overall consistency. Then 
a new value of y was chosen and the process was 



repeated until satisfactory agreement with experiment 
was achieved. No attempt to specify uncertainties 
in the parameters was undertaken. The final choice 
of values for r, y, k, and f is listed in table 2. 



Table 2 



T-Tc 


k, m ! 


T, S 


y 


0.05 
0.1 
1 
2 
5 
10 


2.4 X10 7 
3.71 X 10 7 
1.6 X10 8 
2.48 X 10 8 
4.43 x 10 8 


1.8X10-* 

1.5X10 9 
8X10 10 
6X10" 10 
4x 10 10 


3500 

1630 

129 

72 
32 
12 











The values of the inverse correlation length k, the volume vis- 
cosity relaxation time t and the ratio of the specific heats y as deter- 
mined by fitting eq (12) to the spectrum of critical opalescence. 



In determining y and k we are guided by the expecta- 
tion that these quantities would be proportional to 
some power of T — T c ; namely [10], 



and 



k = A(T-T c )» 

y = B{T-T c )~y +a . 



(20) 

(21) 



The use of y in eq (21) follows standard notation 
for exponents. Although this double meaning for 
y is potentially confusing it should be clear from 
the context whether the exponent y or y= C p /C v is 
intended. The exponents v, y and a specify the tem- 
perature dependence of the correlation length k" 1 , 
the specific heat at constant pressure C p and the 
specific heat at constant volume C v respectively. 
In the interval 0.05 K ^ T— T c ^ 5 K we obtain agree- 
ment with the experimental data using ^ — 0.64, 
^ = 1.6Xl0 8 7&-\ y-a-1.1 and B= 129. 

The correlation length exponent v — 0.64 is con- 
sistent with other estimates [11]. There have been 
no direct determinations of k for C0 2 so it is not 
known whether or not A = 1.6 X 10 8 ra -1 is a reasonable 
choice. An independent determination of k by measur- 
ing the intensity of the scattered light as a function 
of wave vector (scattering angle) is needed. 

The specific heat ratio exponent y — a — 1.1 is 
just barely in agreement with the exponents y and a 
determined from equation of state data [12]. The 
coefficient 5 = 129 is to be compared with an "ex- 
perimental" value of y — 120 for T—T c — 1 K. This 
number was obtained by subtracting out the vibra- 
tional contribution to the specific heat ratio obtained 
from equation of state data. 

Although the Brillouin line-width data may be 
fit using a frequency dependent thermal conductivity, 
the data may also be fit using r, the thermal con- 
ductivity relaxation time, set equal to zero. For 
example a set of values for r were obtained which 
may be represented as 



f = 2Xl0- 11 (7 1 -7V)- 1 s 



(22) 
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in the interval 0.05 K^T-T C ^2 K. This form 
has no significance as for 7*— 7V = 0.05 K,r=5x 10~ 10 s 
and t = yield Brillouin line widths lying within the 
uncertainty of the measured value. Until more precise 
Brillouin line-width data become available it will 
not be possible to decide whether or not a frequency 
dependent thermal conductivity is an appropriate 
feature of a dynamical model of critical point density 
fluctuations. 

The values for the relaxation time r may be approxi- 
mately represented as 



T = 8xl0- 10 |r-7V 



1/3, 



(23) 



The value, — 1/3, of the exponent in eq (23) is an un- 
expected result in that a —1/3 power divergence of 
the volume viscosity has not been suggested by theo- 
retical attempts to determine transport coefficients in 
the critical region [6, 7]. 

Two other sets of parameters were considered. The 
first of these had the specific heat ratio exponent 
y — a=l. While an acceptable fit to the experimental 
data was obtained, this set of parameters was rejected 
because y — a= 1 is not consistent with y and a deter- 
mined from equation of state data. The second set of 
parameters had 7 — a =1.2. It was not possible to 
fit the intensity ratio data with the second set. 




1 1 1 



J 1 1 1 1 1 1 i9 



10 



T-T r ,K 



Figure 2. 



Comparison of experimental and computed intensity 
ratios. 

The solid circles are computed values: the error bars indicate the spread of the experi- 
mental points in figure 3 of reference [3]. The solid line, \80\T—T r \- ]tH , was drawn through 
the computed points "by eye". It is, on this plot, indistinguishable from 185\T—T r \ '■"- 
i the least squares fit to the experimental points. 
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FIGURE 1. Comparison of experimental and computed Brillouin 
line widths. 

The circles are the computed line widths. The solid lines indicate the uncertainty in 
the experimental results. This includes both the quoted uncertainty and the spread of the 
data in figure 3 of reference [3|. 

The type of fit obtained for the Brillouin line width 
is shown in figure 1 and for the intensity ratio as shown 
in figure 2. The error band in figure 1 is taken from 
reference [3] and represents experimental uncertainty 
and spread in the data. No uncertainty for the intensity 
ratio is quoted in reference [3], but the scatter in the 
data is such that our computed values are compatable 
with the experimental values. The central component 
line widths are not indicated as the computed values 
were required to agree with the experimental widths. 

4. Discussion 

It is not reasonable to expect the single relaxation 
time model for the volume viscosity to provide an 



adequate description of density fluctuations in the 
critical region. More probably a model with a "distri- 
bution of relaxation times" will be found when im- 
proved data become available, to yield a better 
description [13]. This would require a modification of 
eq (5) with corresponding changes in eqs (9), (10), 
(13), (14), (15), and (16). 

A fairly good test of the single relaxation time model 
can be made by measuring the Brillouin sound speed 
as a function of wave vector. With a single relaxation 
time the dispersion in vb, the Brillouin speed, occurs 
over a fairly narrow range in k. With a distribution of 
times the range in k over which dispersion in vb is 
exhibited increases and the increase in vb with in- 
creasing k is less abrupt. We have determined vb as 
a function of k using the parameters in tables 1 and 
2. These speeds are plotted in figure 3 as functions of 
T—T c for four different values of the wave vector k. 
The limiting curves for C^ and C are also shown. 
The minimum in vb is due to the decrease in Co as 
T—> T c . The wave vector dependence is shown in 
figure 4. There (vb/Co) 2 for different T—T c is plotted 
as a function of k. If a "distribution of relaxation times" 
model were used, the slope of the curves in figure 4 
would be less steep. 

A dynamical model for long wavelength density 
fluctuations in a fluid near the critical point has been 
described by the linearized hydrodynamic equations 
which were modified to allow for nonlocal pressure- 
density relationship and to allow for frequency depend- 
ence in the volume viscosity. The parameters of this 
model may be chosen so that the calculated spectrum 
is in agreement with the observed spectrum. The 
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FIGURE 3. The Brillouin speed v B versus T — T c for several values 
of wavevector. 

Curve 1, C«; curve 2, it = 2.18 X lO^m" 1 (0=166°) curve 3, k= 1.54X 10 7 m-' (0 = 90°); 
curve 4, k= 1 X lO'in" 1 (0 = 54°): curve 5, Ar = 5xlO»nr I (0 = 26°); curve 6,C . 



these quantities enters in a very simple way. Never- 
theless the intensity provides an important overall 
consistency check on the parameters of the model. 

The second comment is that the k 4 term in the 
central component line width is not simply a measure 
of the two particle correlation length. This is seen 
clearly in eq (19). 

A more careful evaluation of the model could be 
made if two sets of measurements were to be per- 
formed. The correlation length, k -1 , has been used as 
an adjustable parameter. It should be determined from 
other light scattering measurements [14]. 

The values of k used in this paper imply a consider- 
able angular dissymmetry in the intensity of the scat- 
tered light for T— TV < 0.1 K. This is just the region 
where multiple scattering becomes important so it is 
not possible to say if our values of k are consistent with 
the measurements discussed in reference [14]. The 
uncertainty in the Brillouin line width measurements 
is too large to permit any meaningful evaluation of the 
frequency dependence of the thermal conductivity. 
The line width measurements should be made with 
greater precision. 




K, I0 7 m~ 



Figure 4. The dispersion in the Brillouin speed as a function of 
wavevector is indicated as (v B /C ) 2 versus k for several values of 
T— T 

Curve l,'r-|f = 0.05 K: curve 2, T-T C = QA K; curve 3, T-T c = l K; curve 4, T~T r 
= 2 K; curve 5, T~T C =5 K. The limiting values {CJCoT 
for curves 3, 4, and 5. 



are shown as horizontal lines 



parameters are consistent with other measurements 
although this could change as more precise measure- 
ments are made. 

Two comments are in order. The first is that near 
the critical point, the intensity ratio is not a measure 
of the ratio of the specific heats. As shown in figure 2, 
calculated values of the intensity ratio are described by 

^^lsoir-rd- 1 - 04 . 

21 H 

To obtain this set of values we used 

y=i29|r-r c |- 1 - 1 . 

Obviously I c \2l H ^ y— 1. The intensity ratio measures 
not only y but also {v R kr) 2 . Unfortunately neither of 
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